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ABSTRACT 

Heavy squark decays into top and charginos or neutralinos could be an unex- 
pected source of top quarks at hadron colliders. A detailed treatment of these 
processes is necessary for a reliable calculation of both the top quark produc- 
tion cross-section and the standard top quark branching ratio in the MSSM. 
Along this line we compute the electroweak corrections to the sbottom decays 
ba — * Xa t within the Yukawa coupling approximation. The calculation of 
these higher order contributions requires renormalization of both the bottom 
squark mixing angle and of tan/?. This type of corrections gives the leading 
order electroweak effects at low and high tan f3. 



A fundamental missing link between the highly successful Standard Model of the 
Elementary Particle Physics and the highly developed status of the experiments is the 
final elucidation of the nature of the spontaneous symmetry-breaking mechanism (SSB). 
Whether the understanding of the SSB will eventually point towards the existence of 
elementary spinless particles (Higgs bosons) or towards the finding of dynamically bound 
fermion-antifermion spinless states, is far from clear at present. Ultimately, experiment 
must decide. In the meanwhile, the Minimal Supersymmetric extension of the Standard 
Model (MSSM) remains immaculately consistent with all known high precision experi- 
ments at a level comparable to the SM This fact alone, if we bare in mind the vast 
amount of high precision data available at present both from low-energy and high-energy 
physics, should be considered at least as greatly notorious. To this fact we have to add 
another prominent feature, to wit: that the MSSM offers a starting point for a successful 
Grand Unified framework where a radiatively stable low-energy Higgs sector can survive. 
Putting things together it is well justified, we believe, to keep alive all efforts on all fronts 
trying to discover a supersymmetric particle. Undoubtedly, the next Tevatron run, and 
of course also the advent of the LHC, should offer us a gold-plated scenario for finding a 
hint of SUSY, if it is there at all. 

Sparticles not much heavier than a few hundred GeV could be produced in significant 
numbers already at the Tevatron. For instance, selectron production was advocated in 
Ref. to explain a purported non-SM event in the Collider Detector at Fermilab (CDF). 
Subsequently, in Ref. [[| it was argued that half of the top quarks at the Tevatron might 
come from gluino decays into top and stop, g — > tt\. Similarly, we may envision the 
possibility that sbottom squarks are pair produced by the usual Drell-Yan mechanism 
and then decay into top quark and charginos: b a — > xT t- Indeed, this would be the lead- 
ing two-body decay if gluinos are heavy enough that the strong decay mode b a — > g b is 
kinematically blocked upQ. The observed cross-section is then equal to the Drell-Yan pro- 
duction cross-section convoluted over the parton distributions times the squared branching 
ratio. However, in the framework of the MSSM, we rather expect a generalization of this 
formula in the following way (schematically), 



We assume that gluinos are much heavier than squarks, so that their contribution to this 

cross-section through qq — > gg followed by g — > tt\ is negligible. From eq.(P we see 

that, if there are alternative (non-SM) sources of top quarks subsequently decaying into 

1 Squark decays have been discussed at the tree-level in several places of the literature. See e.g. Ref.[Q 
for some relatively recent references on the subject. 
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the SM final state, W + b, then one cannot rigorously place any stringent lower bound on 
BR(t — > W + b) in the MSSM from the present FNAL data. Indeed, we could as well 
have non-SM top quark decay modes, such as e.g. t — > t a Xa an d t — > ||, 
that could serve, pictorially, as a "sinkhole" to compensate (at least in part) for the 
unseen source of extra top quarks produced at the Tevatron from sbottom pair production 
(Cf. eq.(^)). However, one usually assumes that BR{t — > W + b) > 50% in order to 
guarantee the (purportedly) standard top quark events observed at the Tevatron. Thus, 
from these considerations it is not excluded that the non-SM branching ratio of the top 
quark, BR(t — >"new"), could be as big as the SM one, i.e. ~ 50%. We stress that 
at present one cannot exclude eq. ([!]) since the observed t — > W + b final state involves 
missing energy, as it is also the case for the decays comprising supersymmetric particles. 
In particular, if tan/3 is large and there exists a relatively light chargino with a non- 
negligible higgsino component, the alternative mechanism suggested in eq.(|]) could be a 
rather efficient non-SM source of top quarks that could compensate for the depletion in 
the SM branching ratio. 

While the squark production cross-section has already received some attention in the 
literature at the level of NLO radiative corrections ]7J , an accurate treatment of the decay 
mechanisms is also very important to provide a solid basis for experimental analyses of 
the top quark production in the MSSM. Thus in this paper we consider the computation 
of the leading supersymmetric electroweak (SUSY-EW) quantum effects on b a — > Xi t, 
namely the ones induced by potentially large Yukawa-couplings from the top and bottom 
quarks (normalized with respect to the SU(2)l gauge coupling): 

A = — - mt X = — - mb (2) 

''g'V^Mwsmp ' b ~ g ~ V2M W cos/3 ' 

At large (> 20) or small (< 1) tan (3 || these effects could be competitive with the known 
QCD corrections . Since in these conditions the full MSSM quantum effects can be 
rather large, their calculation is indispensable to account for the observed top quark 
production cross-section ([!]) in the MSSM or, alternatively, to better assess how much the 
determination of the SM branching ratio BR(t —>■ W + b) is affected in the MSSM context 
after plugging in the experimental number on the LHS of eq.(|l]). 

We address the calculation of the one-loop corrections to the partial width of b a — > Xi t 
within the context of the on-shell renormalization framework [|lOfl . However, in SUSY 
extensions of the SM, there are additional ingredients concerning the renormalization 
program that must be fixed. Consider the MSSM interaction Lagrangian involving the 
top-sbottom-chargino vertex in the mass-eigenstate basis: 

C th = -g Kxt(A a ;P L + t, t A m P R )t + h.c. , (3) 
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where is the sign of the ith chargino eigenvalue Mj (i = 1, 2 with \M\\ < \M%\) in the 
real matrix representation, and the coupling matrices are denoted byf] 

A% = R la V a - \ b R 2a V i2 , A a 2 = -R la X t U l2 . (4) 

The explicit appearance of the Yukawa couplings @ in the above Lagrangian requires 
both the introduction of top and bottom quark mass counterterms (in the on-shell scheme) 
and also a suitable prescription for tan/3 renormalization. We denote by mi (a = 1,2), 
with mi < mg , the two sbottom mass eigenvalues. The sbottom mixing angle 6~ b is 
defined by the transformation relating the weak-interaction (b' a = 6^,6^) and the mass 
eigenstate (b a = 61,62) squark bases: 

K = Ra b k; R=( C °^ - 8i "M , (5) 
a ' y sm^ cos9~ h ) ' w 

where R is the matrix appearing in eq.(§). From this change of basis the sbottom mass 
matrix, 

M 2 = ( M? l + m 2 b + cos2^{-l + I s 2 w ) M 2 Z m b {A b - // tan (5) \ 

b \ m b {A b - fit&nfi) Ml R +m 2 h -\ cos2^s^M|, ) ' W 

becomes diagonalized as follows: R^ Ai? R = diag jm~ 2 , m^j . 

Formally, we start the renormalization procedure as usual |I| i.e. by introducing 
parameter and field renormalization constants: 

9 

. I TOO 

9 > 



m, — > (1 + - — -) m 

m 



t - (1 + -^^ + -^^)*, 

6 a - (l + ^)6 a + 5rt (6^ a), 

_ 5 tan/? 

tanp — > (1 H — Jtanp, 

tanp 

R -> R + 5R. (7) 

Notice that the sbottom wave-function renormalization involves the mixed scalar field 
renormalization constant (6 7^ a) 

S afe (m? ) 

5Z afc = , \ , (8) 
m; — m~ 



2 See Refs.p, O for full notation niceties. 
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where S afc is the bare mixed self-energy connecting the physical states b a and bj, (a ^ b). 
The renormalization of the Yukawa couplings (0) involves the mass counterterms and the 
counterterm to tan f3 through the relations 

Xt 



5m t 



5M W tan 3 

+ sin p 



M w 
SM W 



— cos 2 P 



tan (3 
8 tan p 



(9) 



m t Mw ' tan/? 

Furthermore, as the last row of eq.(0) shows, we have introduced a counterterm SR for 
the rotation matrix R that diagonalizes the sbottom mass matrix. This counterterm can 
be replaced (see below) by a counterterm 56^ for the sbottom mixing angle. Indeed, we 
shall proceed this way as we wish to take 0g as an input in our renormalization program. 
Substituting the relations (0) into the bare Lagrangian ([3]) we find: 



£onc-ioop = -9 Kxt [A a ; + 5C a + l ) P L + t t [A m + 5C?) P R 



t + h.c. 



(10) 



with 



6Ci 

SC ai 



A a l ( — + -5Z a + -5Z R + -5Z\ \ + 5A% + 5Z ba A b ; , 
\ 9 2 2 2 J 

A™ ( 5 JL + -5Z a + -5Z\ + -8Z R ) + 5A m + 5Z ba A b l . 
\ g 2 2 2 1 



'IV 



The full structure of the four on-shell renormalized decay amplitudes for b a — > Xi t (a = 
1, 2; i = 1,2) follows from the previous Lagrangian after including the contributions from 
the (LH and RH) one- loop vertex form factors Flr' 



i T(b a ->tXi) =igut 



e l (A° 



(12) 



where 



A , 



8C™ + €i F\ 



L i 



Let now Tq* be the tree-level partial width of the decay b a 
Lagrangian (H): 

\(a,i,t) {[(A a ;) 2 + (A a S 



A* = 5Cf + F% . (13) 
Xl t obtained from the 



16 7rm? 

ba 



[mi 



Mf 



AAfA^mteiMi} , 



(14) 

with X(a,i,t) = A(m~ , M?,m|) the usual Kallen function for the given arguments. The 
quantum correction to V™ can be described in terms of the quantities S ai = (T ai — Tg 2 ) /Tq*, 
where T ai is the corresponding one-loop corrected width. From the previous formulae, 5 at 
can be worked out as follows: 



2 (mf - - m 2 ) (A* Af + A^! Ag) - 4 m t e 4 Afj (A*A£ + A™ Aj 
[{Af) 2 + (A™) 2 } (mf - M 2 - m 2 ) - 4 Ae.* m t M< 



(15) 
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For the further procedure, we have to fix the renormalization conditions that determine 
the counterterms £tan/3 and S9 b which enter the explicit structure of S A± in eg . (|Tl~D . 
Following Ref.p we use the decay process H + — > r + u T to define the parameter tan/? 
through 

r <* + - ^> - ^ £ s$ £ tm " 0ii - ArMSSM) ■ (i6) 

This determines the counterterm 5 tan (3 § as follows: 
5 tan 1 /5M^ 5g 2 \ 1 



tan/3 2 1 Mjl g 2 J 2 



-5Z H ± + cot (3 5Z HW + A r , (17) 
l w y~ I 1 

where A r stands for the full set of MSSM one-loop corrections to the partial width of 
H + — > r + v T . For fixing 56^, we require that the renormalized sbottom mixing angle (that 
we use as an input data, Cf. eq.(pO[) below) does not feel a shift from the mixed sbottom 
bare self-energies S afe between the physical states b a and b^ (a ^ b). This is similar to 
the prescription adopted in Refs.|| jl2"| , though it is not identical. In our formalism, the 
3-point Green functions explicitly incorporate the mixed field renormalization constants 
SZ ab (a b) and are therefore renormalized also in the Q~ h parameter. The UV-divergent 
parts of the 3-point functions are cancelled against 89~ b by defining the latter as follows: 



1 iS> )+E 12 (m ) 

86~ b = \ {6Z» - 62P) = \ —^4 A-^ . (18) 

2 2 mf — mi 

Of course another equivalent choice could just be 56~ b = 5Z 12 (or —5Z 21 ), but eq. (p~8"D , is 
more symmetrical; the numerical differences among the finite parts of the two choices are 
negligible. 

The analytical formulation developed thus far is well suited to tackle the general 
problem of the SUSY-EW corrections to squark decays. Since the dominant part is from 
the Yukawa sector we wish to pursue our calculation in the following within the Yukawa 
coupling approximation. This means that we are going to compute the leading electroweak 
effects of 0{Xf) and 0(X 2 b ) that emerge for large values of the Yukawa couplings (0) when 
the remaining gauge contributions - of 0(g 2 )- are subdominant. In practice we shall 
only explore the large tan/3 regime, typically tan/3 > 20; the possibility tan/3 < 1 is 
not so appealing from the theoretical point of view. Thus within our approximation 
we will include the correction A T in leading order 0(X 2 ) of the r Yukawa-coupling, A T . 
Notice furthermore that for Af, >> 1 the tree- level decay rate, eq. (Q), is maximized. 
Therefore, the large tan/3 range is expected to be the most relevant one for the decay 
under consideration. 

In our approach, we set the SU{2)l gaugino mass parameter M » |/x|,Mw in the 
chargino mass matrix [EL and therefore the chargino xt is mainly higgsino, whereas the 
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chargino xf is mainly gaugino and does not contribute to our decays. It is only in this case 
that the Yukawa-coupling approximation makes sense. Thus, since m t - > 80 — 90 GeV, 
in this approach the decay into stop and neutralino t —>■ t a \° a is kinematically forbiddenP]. 
Therefore, in this approximation the relevant counterterms in eq.(p]) boil down to 

8A a + = -5R 2a h - Ria5\b 

6 A"? = -5R la X t -R la 5X tl (19) 

with 5\b and 5X t as in eq.(||). The one-loop Feynman diagrams contributing to b a — > \i t 
in our Yukawa-coupling approximation are depicted in Fig. 1. Specifically, in Fig. la we 
show the wave-function renormalization graphs from which the field renormalization con- 
stants 5Z a , 5Z l and SZ t LR are straightforwardly computed; and in Fig. lb we display 
the contributions going to the vertex form factors F™ H . We refrain from writing out 



cumbersome formulas in this letter [13| . Nevertheless we remark that the calculation is 
well-defined in the Yukawa coupling approximation. In this respect, we have checked the 
effective cancellation of divergences among the bare diagrams and the various counter- 
terms. From the practical point of view we emphasize that the definition of tan f3 that we 
have used to fix 8 tan (3 (Cf. eq.(|l(J)) is quite convenient as the Tevatron (and the LHC) 
collaborations have started to develop all the necessary techniques for r-tagging and if ± 



identification [14 



Let us now pass on to the numerical analysis. It is summarized in Figs. 2-5. The set 
of independent parameters that we use in the sbottom sector consists of the masses and 
the mixing angle 

(m biJ m hj e- b ), (20) 
whereas for the stop sector we just have in addition 

(m~ tv 9i) (21) 

since by SU(2)i gauge invariance the value of the other stop mass mi is already deter- 
mined. Similarly, the sbottom and stop trilinear terms Ab and A t are fixed by the previous 
parameters as follows: 

9 9 9 9 

YYl% — Tf\~ TYl~ — TYl~ 

A b = tan (3 + ^ ^sin2%; A t = ucot/3 + — ^ ^sin2%. (22) 

2 m b 2m t 

We impose the approximate (necessary) condition 

A\ < 3 (mf + m\ + M 2 H + /i 2 ) , (23) 

where m^ is of the order of the average squark masses for q = t,b, to avoid colour-breaking 
minima in the MSSM Higgs potential ||15|| . 



3 See, however, Ref.|L3| for a more general discussion. 
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For the numerical analysis, we assume that the sbottom masses are not too heavy 
( m bi — 350 GeV) and that the charged Higgs mass, M H , is such that the decay t — > H + b 
is available. In this way the quantum corrections that we shall evaluate could have an 
impact already for the Tevatron physics of top quark production. However, in Ref . [[jj 



we include the corresponding study of the radiative corrections to sbottom decays also 
in the case of very heavy ( "obese" ) sbottom and charged Higgs bosons (rriq > 400 GeV, 
Mn > m t) which should be relevant for LHC Of course, tan/? and the SUSY Higgs 
mixing parameter \i are also additional independent inputs for our calculation. In the 
relevant large tan ft segment under consideration, namely 

20 < tan/? < 40, (24) 

the bottom quark Yukawa coupling A& is comparable to the top quark Yukawa coupling, 
X t . Even though the extreme interval 40 < tan/3 < 60 can be tolerated by perturbation 
theory, we shall confine ourselves to the moderate range (pip . This is necessary to preserve 
the condition (|23|) for the typical set of sparticle masses used in our analysis. In particular, 
notice that with our chosen set of independent inputs in eqs.(E0)-(ET|) the A-parameters 
are not free -they are fixed by eq. (|22|) . We point out that the colour stability requirement 
([23]) could be satisfied independently of tan ft if the A-parameters would be chosen directly 
as a part of the set of inputs and then taken sufficiently small. Nevertheless this possibility 
is not so convenient in our analysis where the sparticle masses are the natural inputs that 
we wish to control in order to make sure that sparticles can be produced and decay at 
the Tevatron as explained in connection to eq. ([!]). 
In Figs. 2a and 2b we evaluate the branching ratio 

BK (b a -» Xi t) - ^TjT 

EiTi = S Q r (6 a -> x° a b) + T (b a -> Xi t) 

+ Toib^H-i^ + ^Toiba^^k), (25) 

(where $° = h°, H°, A ) as a function of tan/3 and mg , respectively, in order to illustrate 
the relevant intervals for these parameters where to compute radiative corrections to our 
decay. In these figures we have also indicated the ranges of tan ft and mg excluded by the 
condition (|23|) for the given values of the other parameters. The corresponding corrections 
5 ai fljjD are shown in Figs. 3a and 3b as a function of the lightest stop and sbottom masses, 
respectively. The allowed range for the sbottom and stop mixing angles is conditioned by 
the upper bound on the trilinear couplings and is obtained from eqs.(p2|) and (0). In the 
case of 8~ b , it is rather narrow, so that the physical sbottom masses basically coincide with 
the LH and RH electroweak eigenstates. In the physical 9~ b range, the variation of the 
correction (|I3|) is shown in Fig. 4a. On the other hand the permitted range for the stop 
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mixing angle, is much larger and we have plotted the corrections within the allowed 
region in Fig. 4b. Notice that the sign of the quantum effects changes within the domain 
of variation of 9%. Finally, we display the evolution of the SUSY-EW effects as a function 
of tan (3 (Fig. 5a) and of \x (Fig. 5b) within the region of compatibility with the constraint 
(Effi)- We remark that for \\x\ > 120 GeV and tan j3 > 20 the corrections can be above 20% 
near the phase space limit of the lightest sbottom decay. However, this effect has nothing 
to do with the phase space exhaustion, which is described by the kinematic function 
X(a,i,t) on the RHS of the tree-level expression fll4|), but rather with the presence of 
the dynamical factor in brackets on that equation which also goes to the denominator of 
5 in eq.(|l5T). That factor is fixed by the structure of the interaction Lagrangian of the 
sbottom decay into charginos and top; and, for the parameters in Fig. 5, it turns out to 
vanish near (actually past) the phase space limit in the case of the lightest sbottom (bi) 
decay. However, this is not so for the heaviest sbottom (62) as it is patent in the same 
figure. 

A few more words are in order to explain the origin of the leading electroweak effects. 
One could expect that they come from the well-known large tan j3 enhancement stemming 
from the chargino-stop corrections to the bottom massQ (Cf. Fig. la). Nonetheless this 
is only partially true, for in the present case the remaining contributions in Fig. 1 can be 
sizeable enough. To be more precise, in the region of the parameter space that we have 
dwelled upon the bottom mass contribution is seen to be dominant only for the lightest 



sbottom decay and for the lowest values of tan/? in the range (|24T ). This is indeed the 
case in Fig. 4b where tan (3 = 20 and therefore the bottom mass effect modulates the 
electroweak correction in this process and 5 becomes essentially an odd function of the 
stop mixing angle. On the other hand, from Fig. 5a it is obvious that the (approximate) 
linear behaviour on tan f3 expected from bottom mass renormalization becomes completely 
distorted by the rest of the contributions, especially in the high tan/? end. In short, the 
final electroweak correction cannot be simply ascribed to a single renormalization source 
but to the full Yukawa-coupling combined yield. 

In general the SUSY-EW corrections to T(b a — ► xi t) are smaller than the QCD 
corrections . The reason why the electroweak corrections are smaller is in part due to 
the condition (|23| ) restricting our analysis within the tan/3 interval (|24T). From Figs. 4 
and 5a it is clear that outside this interval the SUSY-EW contributions could be much 
higher and with the same or opposite sign as the QCD effects depending on the choice of 
the sign of the mixing angles^. Moreover, since we have focused our analysis to sbottom 
4 See e.g. || and references therein. 

5 In the conditions of Ref . || the evolution of the QCD corrections with tan f3 would be very slow 
because the assumed sbottom mixing angle is zero. In the present case, however, Q~ h ^ and the QCD 



corrections evolve significantly with tan/3 but become saturated 1 13 
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masses accessible to Tevatron, again the theoretical bound (^3|) severely restricts the 
maximum value of the trilinear couplings and this prevents the electroweak corrections 
from being larger. This cannot be cured by assuming larger values of Mh and/or of /i 
due to our assumption that t — ► H + b is operative and because \i directly controls the 
value of the (higgsino-like) chargino final state in our decay, so that basically we have 
\fi\ < m~ ba —Mn. The restriction cannot be circumvented either if we assume larger values 
of mi a , for it has been shown that too heavy stops are incompatible with the CLEO data 
on b — > S7 both at low and high tan /3 [|I6| , [OJ . We point out that the MSSM analysis 
of b — > s 7 also motivated the sign choice A \x < in our numerical calculation |17| . 
Admittedly, the situation with radiative S-decays is still under study and there are many 
sources of uncertainties that deserve further experimental consideration. Still, we have 
used this information to focus on a limited domain of the MSSM parameter space. 

Finally, we mention that the treatment of the electroweak corrections to the closely 
related decay t a — > %° t is similar to the one presented here, but for lack of space it will 
be presented elsewhere |TJ| . It suffices to say that the branching ratio and the radiative 
corrections to that decay can be sizeable only for small values of tan /3, so that this 
process becomes optimized in a region of parameter space complementary to the one 
under consideration. 

In summary, the MSSM corrections to squark decays into charginos and neutralinos 
can be significant and therefore must be included in any reliable analysis of top quark 
physics at the Tevatron within the MSSM. The main corrections stem from the strongly 
interacting sector of the theory (i.e. the one involving gluons and gluinos), but also 
non-negligible effects may appear from the electroweak sector (characterized by chargino- 
neutralino exchange) at large (or very small) values of tan/?. Failure of including these 
corrections in future studies of top quark physics at the Tevatron, both in the production 
and decay mechanisms, might seriously hamper the possibility of discovering clear-cut 
traces of SUSY physics from the identification of large non-SM quantum corrections in 
these processes. As already stated, we have mainly concentrated on the impact of these 
quantum signatures in the physics of the Tevatron, but important effects are also ex- 



pected [|13|] for experiments aiming at the production and decay of "obese" squarks at the 
LHC The latter type of squarks could be free of some of the restrictions that have been 
considered for the present calculation. 
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Figure Captions 

• Fig.l SUSY-EW Feynman diagrams, at one-loop order, correcting the partial width 
of b a — > Xi t'- ( a ) Self-energy diagrams; (b) Vertex contributions. There are other 
possible diagrams, but they do not contribute in the Yukawa-coupling approxima- 
tion. 

• Fig. 2 (a) The branching ratio of b a — > Xi t as a function of tan/? for the various 
decays a = 1, 2 with < m^; (b) As in (a), but as a function of . The marked 
parts of the abscissa in both figures are excluded by the condition (|23|). The fixed 
parameters for (a) and (b) are given in the frame. 

• Fig. 3 (a) The SUSY-EW corrections (|T5|) as a function of m^; (b) As in (a), but 
as a function of m~ bi . Rest of inputs as in Fig. 2. 

• Fig. 4 (a) Evolution of the SUSY-EW corrections as a function of the sbottom 
mixing angle, 8^, within its allowed range; (b) As in (a), but as a function of the 
stop mixing angle, % Remaining inputs are as in Fig.2 

• Fig. 5 (a) The SUSY-EW correction as a function of tan/3; (b) As in (a), but as a 
function of //. Rest of inputs and notation as in Fig.2. 
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